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trigonal and hexagonal substrates, in which misfit disloca- 
tions are reduced by approximate lattice matching of the 
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substrate structure, enabling the development of alloyed 
semiconductor layers of greater thickness, resulting in a new 
class of semiconductor materials and corresponding devices, 
including improved hetero -bipolar and high-electron mobil- 
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Atomic alignment which was interpreted from 
XRD data shows 30 degree rotation: Type [B] 



Sapphire, view on {0001} 




Sapphire (AI2O3) 
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Fig. 2 
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Fig. 3 



Coordinates are rotated by 30° but the angle between 
reciprocal basis vectors, Clmti 3 is only 60° 
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Fig. 4a Fig. 4b Fig. 4c 
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Fig. 5 
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EPITAXIAL GROWTH OF CUBIC 
CRYSTALLINE SEMICONDUCTOR ALLOYS 
ON BASAL PLANE OF TRIGONAL OR 
HEXAGONAL CRYSTAL 

CROSS-REFERENCE TO RELATED 
APPLICATIONS 

This application claims the benefit of the respective filing 
dates of the following commonly assigned U.S. Provisional 
Patent Application Ser. No. 60/980,870 filed on Oct. 18, 2007, 
Ser. No. 60/980,871 filed on Oct. 18, 2007, Ser. No. 60/980, 
878 filed on Oct. 18, 2007, Ser. No. 60/980,880 filed on Oct. 
18, 2007, and Ser. No. 60/980,881 filed on Oct. 18, 2007, each 
of which contains an overlap of inventive entity with the 
present application. In addition, this application incorporates 
by reference the entire disclosures of the following com- 
monly assigned nonprovisional U.S. patent applications 
being filed on the same date as the present application Ser. No . 
12/254,016, entitled “THERMOELECTRIC MATERIALS 
AND DEVICES;” Ser. No. 12/254,134, entitled “HYBRID 
BANDGAP ENGINEERING FOR SUPER-HETERO-EPI- 
TAXIAL SEMICONDUCTOR MATERIALS, AND PROD- 
UCTS THEREOF;” Ser. No. 12/288,379, entitled “RHOM- 
BOHEDRAL CUBIC SEMICONDUCTOR MATERIALS 
ON TRIGONAL SUBSTRATE WITH SINGLE CRYSTAL 
PROPERTIES AND DEVICES BASED ON SUCH MATE- 
RIALS;” Ser. No. 12/288,380, entitled “X-RAY DIFFRAC- 
TION WAFER MAPPING METHOD FOR RHOMBOHE- 
DRAL SUPER-HETERO-EPITAXY;” and Ser. No. 12/254, 
150, entitled “METHOD OF GENERATING X-RAY 
DIFFRACTION DATA FOR INTEGRAL DETECTION OF 
TWIN DEFECTS IN SUPER-HETERO-EPITAXIAL 
MATERIALS;” each one claiming priority to the above cited 
provisional applications. 

ORIGIN OF THE INVENTION 

This invention was made in part by employees of the 
United States Government and may be manufactured and 
used by or for the Government of the United States of 
America for governmental purposes without the payment of 
any royalties thereon or therefor. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention is in the fields of solid-state physics, 
semiconductor material and heteroepitaxy, and more particu- 
larly concerns lattice-matching methods for rhombohedral 
growth of crystals comprised of semiconductor alloys of 
group IV and other cubic crystalline materials on the basal 
plane of trigonal or hexagonal crystalline substrates, and the 
resulting semiconductor compositions. 

2. Description of the Related Art 

The alloys of group IV elements, including carbon (C), 
silicon (Si), germanium (Ge), and tin (Sn) are important semi- 
conductor materials. For example, SiGe alloy is widely used 
in Hetero Bipolar Transistors (HBT) and High Electron 
Mobility Transistors (HEMT) for high speed computing and 
wireless communications, as well as in Thermoelectric (TE) 
materials. The growth of these alloys has been very difficult 
because the alloy has a different lattice constant from that of 
the commonly used substrate, such as Si wafer. Prior efforts 
have tried to match the lattice constant of cubic structures 
directly, such as by the approach wherein L target =L substrate 
(the lattice constant of cubic substrate). However, even a 


2 

small mismatch of these parameters can result in extremely 
high pressures in the alloy layer. Thus, in prior efforts, the 
thickness of the grown layers was very small and often had a 
high density of defects such as misfit dislocations. 

5 Similar considerations apply with respect to other cubic 
crystalline semiconductor materials, including group III-V 
and II-VI materials in zinc-blende structure, as well as other 
materials in Body Centered Cubic (BCC) and Face Centered 
Cubic (FCC) structures, 
to 

SUMMARY OF THE INVENTION 

It is an object of the present invention to develop new 
methods for controlling the growth of rhombohedral crystals 
15 composed of group IV alloys and other cubic crystalline 
materials on trigonal and/or hexagonal crystal substrates, and 
thereby to provide new and useful alloyed semiconductor 
materials which previously could not be reliably produced. 

Accordingly, the present invention provides heteroepi- 
20 taxial semiconductor materials comprising rhombohedral 
crystals of group IV semiconductor alloys and other cubic 
crystalline materials grown on the basal plane of trigonal and 
hexagonal substrates in which misfit dislocations are reduced 
by lattice matching of the overlying rhombohedral to under- 
25 lying trigonal or hexagonal substrate structure. This enables 
the development of alloyed group IV and other cubic crystal- 
line semiconductor layers of greater thickness, resulting in a 
new class of semiconductor materials and corresponding 
devices, including improved heterobipolar and high-electron 
30 mobility transistors, and high-mobility thermoelectric 
devices. 

To achieve these results, the present invention generally 
uses a new lattice matching technique to calculate alloy com- 
positions in order to approximately match substrate lattices, 
35 and thereby allow the growth of group IV alloy and other 
cubic crystal layers that are much more free of dislocation 
defects than those that could be obtained under prior 
approaches. The present approach relates to any alloys of 
silicon (Si), germanium (Ge), carbon (C) and tin (Sn) (gen- 
40 erally referred to in this specification as “SiGeCSn”) on any 
trigonal/hexagonal substrate, such as sapphire (i.e. single 
crystalline A1 2 0 3 ) (or hexagonal substrates, such as quartz), 
as well as to other cubic crystalline materials, including group 
III-V and II-VI materials in zinc -blende structure, and poten- 
45 tially to other materials in Body Centered Cubic (BCC) and 
Face Centered Cubic (FCC) structures. 

In one embodiment, lattice matching equations, when sat- 
isfied by a particular alloy composition, will result in good 
crystal growth of the alloy on rhombohedral, trigonal, or 
50 hexagonal substrates. Further, techniques for calculating 
Bragg angles and inter-planar angles for twin detection for all 
group IV alloys as well as group III-V and II-VI materials in 
zinc-blende structure, on any trigonal or hexagonal substrate, 
may be applied to a new X-ray diffraction method for twin 
55 detection developed in accompanying patent application Ser. 
No. 12/254,150, entitled “METHOD OF GENERATING 
X-RAY DIFFRACTION DATA FOR INTEGRAL DETEC- 
TION OF TWIN DEFECTS IN SUPER-HETERO-EPI- 
TAXIAL MATERIALS,” and Ser. No. 12/288,380, entitled 
60 “X-RAY DIFFRACTION WAFER MAPPING METHOD 
FOR RHOMBOHEDRAL SUPER-HETERO-EPITAXY,” 
which have both been incorporated herein by reference (spe- 
cific XRD methods are claimed in the separate applications). 

The present invention makes it possible, for example, to 
65 incorporate a high content of germanium and tin atoms into 
silicon without critical thickness limit because it is lattice- 
matched, thereby providing an increase in the built-in electric 
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field for faster SiGeCSn HBT and HEMT devices and more 
efficient SiGeCSn TE devices than previously possible. Other 
aspects and advantages of the present invention will be appar- 
ent from the accompanying drawings, and the detailed 
description that follows. 5 

BRIEF DESCRIPTION OF THE DRAWINGS 

For a more complete understanding of the present inven- 
tion and the advantages thereof, reference is now made to the to 
following description taken in conjunction with the accom- 
panying drawings, wherein like reference numerals represent 
like parts, in which: 

FIG. 1 graphically depicts in FIG. 1 {a) the primitive unit 
cell of trigonal crystal structure of sapphire (A1 2 0 3 ) and FIG. 15 
1 (b) the diamond structure of Si, Ge, C and Sn. 

FIG. 2 shows the crystal lattice constant relationship 
between rhombohedral transformation of diamond cubic 
structure and simplified crystal structure of sapphire. The 
actual crystal structure of sapphire is very complex and not 20 
shown here. 

FIG. 3 shows plots of two kinds of X-ray diffraction data 
which were obtained at different Euler angles (the two sets of 
data were obtained at Phi angles which were apart by 60°). 

FIG. 4 shows the direction of SiGe [220] in FIG. 4 (a) 25 
real-lattice and sapphire [1-104], in FIG. 4 ( b ) reciprocal 
lattice and in FIG. 4 (c) real-space lattice. All arrows are 
projected on the basal plane. 

FIG. 5 shows atomic alignment of SiGe on sapphire, in a 
Type [B] alignment. 30 

FIG. 6 shows 30° rotated atomic alignment of GaN on 
sapphire, in a Type [B] alignment. 

FIG. 7 shows a comparison between Type [A] and Type [B] 
atomic alignment, with different lattice matching conditions. 

FIG. 8 shows fractional “ 2 A of Type [A] alignment.” 35 

FIG. 9 shows fractional “Vi of Type [B] alignment.” 

FIG. 10 illustrates (Top) lattice constants of group IV semi- 
conductor elements, (Lower Left) fast electron mobility in 
SiGe HBT, and (Lower Right) even faster electron mobility in 
SiGeCSn HBT. 40 

DETAILED DESCRIPTION 

The following is a detailed description of certain embodi- 
ments of the invention chosen to provide illustrative examples 45 
of how it may preferably be implemented. The scope of the 
invention is not limited to the specific embodiments 
described, nor is it limited by any specific implementation, 
composition, embodiment or characterization depicted in the 
accompanying tables or drawings or stated or described in the 50 
invention summary or the abstract. In addition, it should be 
noted that this disclosure describes a number of methods that 
each comprise a plurality of steps. Nothing contained in this 
written description should be understood to imply any nec- 
essary order of steps in such methods, other than as specified 55 
by express claim language. 

In the ensuing description, the well-known Miller indices 
notation of lattice planes is used. That is, crystal planes are 
designated by numbers within “( )”, groups of similar planes 
are designated by numbers within “ { }”, direction or length is 60 
designated by numbers within “[ ]”, and groups of similar 
directions are designated by numbers within “< >”. 

Sapphire has a trigonal crystal structure as shown in FIG. 1 
(< a ). The group IV elements, Si, Ge, C, and Sn, have cubic 
diamond structures as shown in FIG. 1 ( b ) (and group III-V 65 
and II -VI crystalline materials in zinc -blende configuration 
have a similar crystal structure). Therefore, the direct com- 


parison of lattice constants does not have a meaning since the 
cubic basal plane, i.e. {001} planes according to Miller index 
classification, cannot be aligned on the trigonal basal plane, 
i.e. the {0001} or “C” plane. However, if we consider a 
rhombohedral transformation of cubic crystals along the 
<11 1> direction, we can see that trigonal {111} planes of 
cubic crystals can match the trigonal {0001 } planes easily, as 
shown in FIG. 2. In this case, a new lattice parameter of cubic 
crystal, shown as a* in FIG. 2 , can match the corresponding 
lattice parameter, a, of the underlying trigonal or hexagonal 
crystalline substrate. Here, the [220] atomic distance a' is 
given by the rhombohedral transformation relation: 


“Type [A]” alignment is used for the above atomic arrange- 
ment at the interface of epitaxial material and substrate mate- 
rial. 

XRD methods are used to confirm the present crystal 
growth methodology and the resulting materials. FIG. 3 
shows two sets of X-ray diffraction data (angle-dependent 
single crystalline XRD) of one SiGe layer on sapphire sub- 
strate which was grown at a certain growth condition. Two 
sets of data were obtained at the same 2 -Theta (detector), 
Omega (goniometer), Chi (tilt) angles, but at different Phi 
angles, which were apart by 60°. In this tilted XRD data, 
(111), (113), (220), and (224) peaks were obtained in 0-20 
scan, which means that poly-crystalline SiGe films were 
formed on sapphire under a certain growth condition. It is 
important to notice that the first and the second order peaks of 
sapphire (1-104) were obtained in the second set of data, but 
there is no sapphire (1-104) peak in the first set of data. The 
lack of (1-104) peak in the first set of data shows the trigonal 
symmetry of sapphire substrate such that 60° rotational sym- 
metry is broken and only 120° rotational symmetry is valid in 
sapphire. The (220) peak in the first set of data shows one 
domain and another (220) peak in the second set of data 
shows the other domain, which was aligned with sapphire’s 
(1-1 04). The coexistence of SiGe’ s (220) peak and sapphire’ s 
(1-104) peak in one single 0-20 scan shows a different new 
atomic alignment scheme, as shown in FIG. 4 and FIG. 5 . 

In a non-cubic crystal system, the reciprocal lattice orien- 
tation can be different from the real space lattice orientation. 
In a trigonal/hexagonal system, the reciprocal basis vectors 
(shown as a and (3 in FIG. 4 ( b )) have 30° rotation from 
real-space lattice vectors (shown as x and y in FIG. 4 (c)) but 
the angle between the reciprocal basis vectors, a and |3 is not 
120°, but only 60° according to the definition equation below. 
On the other hand, the cubic system’s reciprocal basis vectors 
are the same as real space basis vectors. Therefore, extra-care 
must be paid to the crystal plane orientation between the cubic 
system and the trigonal/hexagonal system because of the 
transformation of the basis vectors. The vector transforma- 
tion relation is given in the following equation (See N. W. 
Ashcroft, N. D. Mermin, Solid State Physics, in Chapter 5 
“The Reciprocal Lattice”, Harcourt College Publishers, New 
York, ISBN#0-03 -083 993 -9 (1976)): 


ai -{a 2 xa 3 )’ 


a { ■( fl 2 Xfl 3 )’ 
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-continued 

_ a\ Xfl? 

b 3 = 2n— — ■— — 7— 
flr(fl 2 Xfl 3 ) 

where ~a x , a 2 , and a 3 are the basis vectors of real- space and , 
b 2 , and b 3 are the basis vectors of reciprocal space. 

In FIG. 4, al, and a2 are shown as x and y, bl and b2 are 
shown as a and (3. The comparison of crystal orientations 
between cubic and cubic, or between trigonal/hexagonal and 10 
trigonal/hexagonal systems do not have this problem. By 
contrast, the comparison of crystal orientations between 
cubic and trigonal/hexagonal systems must consider the rota- 
tion and inter- vector-angle change of lattice basis vectors in 
reciprocal space of the trigonal/hexagonal crystal. Because 15 
the line 401 in FIG. 4 (a) and line 403 in FIG. 4 (c) have to be 
aligned in the same direction for two peaks to appear in the 
same 0-20 XRD scan, the SiGe crystal is atomically aligned 
to sapphire’s hexagonal basal plane as shown in FIG. 5, which 
shows 30° rotation. 20 

In FIG. 5, sapphire substrate is shown as hexagon 501 with 
three dots that represent trigonal symmetry, and SiGe epi- 
taxial layer is shown as superimposed hexagon 503 with three 
tiny parallelograms that is aligned with 30° rotation with 
respect to hexagon 501. Six angle-dependent single crystal- 25 
line XRD data sets are shown in the corresponding azimuthal 
orientation in Euler angles. These six XRD data sets are 
categorized in two groups (505a -c and 507 a-c) having similar 
nature respectively, indicating trigonal symmetry. Therefore, 
all the drawings and data in FIG. 5 show that the grown 30 
epitaxy layer is atomically aligned according to a new atomic 
alignment lattice-matching scheme, which is characterized 
by 30° rotation between the projected lattice axis of epitaxial 
layer and substrate layer. This additional atomic alignment 
will be referred to as “Type [B]” alignment. 35 

Some other material systems also show this kind of atomic 
alignment, which is identified by 30° rotation of rhombohe- 
dral epitaxial layer on trigonal (or hexagonal) substrate. For 
example, it has been reported previously that gallium nitride 
(GaN) on sapphire has 30° rotation between GaN epitaxial 40 
layer and sapphire substrate as shown in FIG. 6. (SeeYamagu- 
chi T, Araki T, Santo Y, et al, Journal of Crystal Growth 237: 
993-997 Part 2 (2002), Akasaka T, Kobayashi Y, Ando S, 
Kobayashi N, App. Phys. Lett. 71 (15): 2196-2198 (1997), H. 

M. Manasevit, F. M. Erdmann, and W. I. Simpson, J. Electro- 45 
chem. Soc. 1 18, 1864 (1971), andT. Sasaki and S. Zembutsu, 

J. Appl. Phys. 61, 2533 (1987).) 

Both atomic alignment configurations, Type [A] and Type 
[B], have to be considered in order to expand the lattice 
matching methods for various group IV alloys and group 50 
III-V and group II-VI crystalline materials in zinc-blende 
structure on many different trigonal or hexagonal substrates. 
FIG. 7 shows the lattice matching conditions in Type [A] 
alignment and Type [B] alignment (See also FIG. 2 with 
regard to Type [A]). The lattice matching conditions are dif- 55 
ferent for Type [A] and Type [B]. See equations (13)-(16) 
below.) 

In addition, comparable alignments are possible with simi- 
lar orientations, where the rhombohedral structures are 
smaller in relation to the underlying trigonal/hexagonal struc- 60 
ture than in the cases of the Type [A] and [B] alignments 
discussed above. This is shown in FIGS. 8 and 9. The addi- 
tional lattice-matching conditions for these possible arrange- 
ments must be considered as well. 

In this reduced form of Type [A] alignment, if the rhom- 65 
bohedral dimensions are reduced relative to the underlying 
trigonal/hexagonal dimensions, there is a third lattice-match- 
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ing condition as shown in FIG. 8, in which the coincidence 
lattice constant, a# (reflecting coincidence with the inner 
vertices of the underlying trigonal/hexagonal lattice), is 2 A of 
the trigonal/hexagonal lattice constant a. We refer to this type 
of alignment as “ 2 A of Type [A]” alignment. 

Similarly, in the reduced form of Type [B] alignment, if the 
rhombohedral dimensions are reduced relative to the under- 
lying trigonal/hexagonal dimensions, there is a fourth lattice- 
matching condition as shown in FIG. 9, in which the coinci- 
dence lattice constant, a" (reflecting coincidence with the 
inner vertices of the underlying trigonal/hexagonal lattice), is 
2,V3 times the trigonal/hexagonal lattice constant a. We refer 
to this type of alignment as “V 2 of Type [B]” alignment 

With these four alignments and lattice matching condi- 
tions, a very widely applicable technology for lattice matched 
rhombohedral growth of group IV alloys has been developed, 
as well as group III-V and II-VI crystalline materials in zinc- 
blende structure, on trigonal or hexagonal substrates. The 
detailed technology is described below for Type [A], Type 
[B] , 2 /3 of Type [A] and Y 2 of Type [B] . 

Lattice Matched Rhombohedral Growth of Cubic Crystalline 

Alloys on Trigonal or Hexagonal Substrates 

[1 ] Properties of Pseudo Alloys in Cubic Crystalline Structure 

In the first order approximation where the material prop- 
erty of alloy changes linearly with composition, the lattice 
constant L of an alloy is determined by: 

L{A^ \-^) = L a +( 1 -x)L b for binary compounds, AJ$ { _ x (3) 

L^AJSyC l _ x _ y )=xL A +yL B +(l-x-y)L c for ternary com- 
pounds, AJByC^ (4) 

L(4^ y C z D i . x ^. z )=xL A +yL B +zL c +(l -x-y-z)L D for 

quaternary compounds, A x B y C z D l _ x . y _ z (5) 

When the alloy has a non-linear second-order term due to 
the mixing entropy, the lattice constant of binary alloy can be 
given by: 

L(x)=A+Bx+Cx 2 ( 6 ) 

where C is called the Bowing Parameter and is four times the 
deviation of L(x) from linearity at x=0.5. 

Similar mathematical relations can be made for non-linear 
properties, such as direct-indirect band-gap transition of 
binary, ternary, and quaternary alloys as well. 

In many cases, the lattice constants of most of the group IV 
alloys (as well as group III-V and II-VI crystalline materials 
in zinc -blende structure) can be approximated with simplest 
linear approximation without using Bowing parameters, 
because the Bowing parameters of lattice-constants are usu- 
ally small. On the other hand, the Bowing parameters for 
band-gap energy are significant and cannot be ignored. 

For group IV elements with diamond structures, the lattice 
constants are given as: L c (Carbon) 3.57 A, L Si (Silicon)=5.43 
A, L Ge (Germanium) 5.66 A and (Grey Tin diamond 
structure)=6.48 A. For example, the lattice constant of 
Si 0 2 Ge 0 3 C 0 ^Sn^ can be calculated to be 0.2xL Sz +0.3x 
L c?e +0.1xL c +0.4xL Sw =5.73 A. 

Based on the foregoing, it is possible to change the lattice 
constant of a pseudo -alloy to a desired value in order to match 
it to the corresponding atomic distance of the substrate. How- 
ever, a few percent change in the lattice constant may be 
equivalent to extremely high pressure in the alloy layer. 
Therefore, an unmatched alloy layer is prone to develop crys- 
talline defects, such as serious threading dislocations that can 
destroy electronic devices. 

For binary alloys, lattice-constant matching in accordance 
with the above equations can occur at or about only one 
alloy-composition since it has only one variable (x). On the 
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other hand, ternary alloys have two variables (x and y) and 
quaternary alloys have three variables (x, y, andz). Therefore, 
lattice matching can occur at or about different compositions 
while varying multiple variables. This additional degree of 
freedom can control additional physical properties, such as 
band-gap energy and thermal expansion while the alloy is still 
in the lattice-matching condition. The following equations 
show the method to fix the composition to match multiple 
physical properties, such as index of refraction, band-gap 
energy, lattice constant, absorption coefficient, work-func- 
tion, and so on. 


Therefore, L target should equal the transformed atomic dis- 
tance, such as: 


Ltarget = ~j= aTri/Hex for Type [A] alignment 
V3 


(14) 


As another example, in the type [B] alignment, the follow - 
io ing equation holds as shown in FIG. 4. 


Three properties will now be considered. L is a lattice 
constant, E is band-gap energy, a is a linear thermal expan- 
sion coefficient. 

(1) For Binary Alloy (AJB 1<C ): 

xxL a + ( 1 -x)xL b =L T arg et : Matching the lattice constant (7) 

(2) For Ternary Alloy (A^B^Cj^): 

xxL A +yxL B +(l-x-y)xL c =L Target : Matching the lattice 

constant (8) 


, _ VT 2 (15) 

a [210] - ~ 2 ~ ^target ~ a Tri/Hex ~ a Tri/Hex 

15 

where a* is the atomic distance of [2-20] direction in FIG. 2, 
^ Trig/ Hex i s th e lattice constant of basal plane of the trigonal or 
hexagonal crystalline substrate and Trig/Hex * s the trans- 
20 formed atomic distance that is aligned to a'. 

Therefore, L to , should equal the transformed atomic dis- 
tance, such as: 


xxE A +yxE B +(\ -x-y)xE c =E Target \ Matching the band- 

gap energy (9) 

Two equations with two variables give one answer with 
unique x and unique y. Only the answer with x and y between 
0 and 1 is meaningful. This answer matches two physical 
parameters at the same time. 

(3) For Quaternary Alloy A*B y C£> Ux _y_ z y 

xxL A +yxL b +zxL c +( 1 -x-y-z)xL D =L Target- Matching 

the lattice constant, (10) 

xxE A +yxE B +zxE c +(l-x-y-z)xE D =E Target- Matching 

the band-gap energy, and (1 1) 

x+a A +yxa B +zxa c +(l-x-y-z)xa D =a Targe t' Matching 

the linear thermal expansion coefficient. (12) 

Three equations with three variables give one answer with 
unique x, unique y, and unique z. Only the answer with x, y, 
and z between 0 and 1 is meaningful. This answer matches 
three physical parameters at the same time. 


Ltarget = ^Tri/Hex for Type [5] alignment 

As noted above, two other fractional coincidence atomic 
alignments are possible as shown in FIGS. 8 and 9, i.e., “Two 
30 Thirds ( 2 A) of Type- A” and “Half (V 2 ) of Type-B”. 

For 2 A of Type [A] alignment as shown in FIG. 8, the 
following equation holds: 


, _V3 _ „ _2 d 7 ) 

a l220l ~ 2 ^ Jtar S et ~ a Tri/Hex ~ ^ a Tri/Hex 

where a' is the atomic distance of [2-20] direction in FIG. 2 
40 and a# Trig/Hex = 2 A-a Trlg/Hex for “% of Type-A” matching and 
a Trig/ Hex is the lattice constant of basal plane of the trigonal or 
hexagonal crystalline substrate. 

Therefore, L tarset should equal the transformed atomic dis- 
tance, such as: 

45 


Lattice Matching Technique 

As noted above, previous attempts were made to match the 
lattice constant of cubic structures directly by an approach 
wherein L targBt =L substrMe (the lattice constant of cubic sub- 
strate). In the present approach, L target should aim to match a 
corresponding atomic distance of a trigonal/hexagonal sub- 
strate, not the direct lattice constant of a cubic substrate. To do 
this, equation (14) below for type [A] alignment and equation 
(16) below for type [B] alignment are used. The following is 
the derivation of these equations. 

For example, in the type [A] alignment, the following 
equation holds as shown in FIG. 2. 


, VT, 

^[ 220 ] — 2 J~* ar 8 et ~ a Tri/Hex 


(13) 


4 (18) 

Ltarget = ^j=a Tri/Hex for“2/3 of Type [A]” alignment 


50 


For the “half of type [B]” alignment, the following equa- 
tion holds, as shown in FIG. 4: 


55 , _ V3" ^ _ 1 „ 1 

a [220l - 2 ^t Qrget ~ 2 aTri / Hex ~ a Tri/Hex 


(19) 


where a' is the atomic distance of [2-20] direction in FIG. 2, 
60 a rngiHex 1S lattice constant of basal plane of the trigonal or 

hexagonal crystalline substrate and 


1 „ 

2 a Tri/Hex 

where a' is the atomic distance of [2-20] direction in FIG. 2 65 
and & Trig/ Hex i s ^e lattice constant of basal plane of the trigo- 
nal or hexagonal crystalline substrate. 


is the transformed atomic distance that is aligned to a'. 
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Therefore, L target should equal the transformed atomic dis- 
tance, such as: 


= 3 07 W///CC for “1/2 of Type [5]” alignment 

There are many crystal wafers that have a trigonal crystal 
structure. The following list shows a few selected materials 
with their lattice constants, a Tri and a Hex , as the lattice con- 
stants of basal plane and c as vertical lattice constant. 
Trigonal Crystals (by Point Group Symmetry): 

Sapphire (A 1 2 0 3 ), a rrz =4.758 A, c=12.99 A 
Gallium Nitride (GaN), a rw =3.189 A, c=5.185 A 
Aluminum Nitride (AIN), a rrz =3.1 1 A, c=4.98 A 
Indium Nitride (InN), a rrz =3.54 A, c=8.70 A 
Zinc Oxide (ZnO), a rrz =3.25 A, c=5.21 A 
Lithium Niobate (LiNb 0 3 ), a rrz =5.148 A, c=13.863 A 
Lithium Copper Tantalate (LiCuTa03, Li=0.667, Cu=0.333) 
a r „.=5.2004 A, c=13.788 A 

Lithium Tantalum Nitride Oxide(LTNO) (Lil6N801Ta2), 
a rrz -=5.983 A, c=9.168A 

(In this list, trigonal lattice was described in the hexagonal 
frame instead of the rhombohedral frame.) 

Hexagonal Crystals (by Point Group Symmetry): 

Quartz, a Hex =4.91 A, c=5.40 A 

Boron Nitride (BN) in Hexagonal form, (Space group D 6C - 
P 6 3 mmc), a i7ex =2.50 A, c= 6.66 A 

Therefore, the equations (14), (16), (18) and (20) can be 
applied to all trigonal crystal and hexagonal crystal wafers 
made with the materials, including but not limited to the 
above list. The use of these equations is not limited to group 
IV alloys in diamond structure. They may also be applied to 
group III-V and II-VI materials in zinc-blende structure. The 
above lattice matching equations (14), (16), (18), and (20) 
apply for all of these materials. The trigonal and hexagonal 
crystals can both be described in hexagonal frame, but the 
main difference between the two structures is that the hex- 
agonal crystal allows 60° rotational symmetry, whereas the 
trigonal crystal has only 120 ° rotational symmetry without 
60° rotational symmetry. 

With this general expansion of technology, we formulate a 
group IV alloy made with any of Si, Ge, C, and Sn, or group 
II- V and II-VI materials in zinc-blende structure, on any of 
the trigonal or hexagonal substrates, whose alloy composition 
is designed to hold the equation (14) in Type [A] alignment, 
equation (16) in Type [B] alignment, (18) in 2 A of Type [A] 
alignment, or equation (20) in half of Type [B] alignment 
within ± 1 0% error of alloy composition. We also note that the 
same analysis should be applicable to other materials in Body 
Centered Cubic (BCC) and Face Centered Cubic (FCC) 
structures as well. 

For example, Si 0 14 Sn 0 86 can match lattice constant of 
sapphire in Type [B] alignment. Therefore the alloy of S^Sn^ 
on sapphire in the range of y=0.86±0.086(10%)= 
0.774-0.946 and x=l-y =0.226-0. 05 4 is within this formu- 
lation. 

General XRD Method for Cubic Crystalline Material on any 

Trigonal/Hexagonal Substrate 

(1) Bragg Angle of X-Ray Diffraction Peak: 

Generally, X-ray diffraction occurs whenever the Bragg 
condition is satisfied: 


2^*-sin0=K-X (b 21) 

where d is the atomic spacing distance of plane (h,k,l), 0 is the 
X-ray diffraction angle, Z is the wavelength of the X-ray, and 


10 

n is an integer. Here d(h,k,l) is determined by the crystal 
structure and the lattice constant, L. 

In cubic diamond structure, 


5 




( 22 ) 


L*^( x ,y,z) is defined to represent the lattice constant of a 
pseudo-alloy such as L (A X B 1-X ) L(A X B^C 1 _ X _ > ,) and L(A X B V 
C z D 1 _ x _ y _ z ), as shown in the previous equation (3), (4), and (5). 
Then, 

15 


d h jc,i(x, y, z) = 


Lqlloyix, y, Z) 
y/h 2 +k 2 + 1 2 ' 


(23) 


20 

Therefore, the XRD peak occurs at the angle (x,y,z), 
such that: 


25 


sin 9hjc,i(x + y + z) 


n-X 

2 ■dh'kjix, y, z ) 


n-X-yJh 2 +k 2 +P 

2 ' Lattoy {x, y, Z) 


(24) 


(2) Inter-Planar Angle: 

3° 

For a cubic diamond structure, the inter-planar angle cubic, 
tycubic’ between two planes, (h 1; k 1? IQ and (h 2 , k 2 , 1 2 ), is given 
by: 

35 

/ii /i2 + k{ k 2 + / 1/2 (25) 

COS (pcubic = , . 

y (h\ -\rk\ + 1\) ■ (hi +kl+ l 2 ) 

40 (For general X-ray diffraction analysis, refer to Elements of 
X-ray Diffraction by B. D. Cullity, Addison- Wesley Publish- 
ing Co., 2 nd or higher edition.) 

Similarly, for a trigonal crystal in hexagonal frame and a 

45 hexagonal crystal, the inter-planar angle, § Tri / Hex , between 
two planes, (h l5 k l5 IQ and (h 2 , k 2 , 1 2 ), is given by: 


50 


55 


COS(p T ri/Hex = 


1 3a 2 

hih 2 + k\k 2 + - (h\k 2 +kih 2 ) + 

\ 

( 

Aj+*f + Ai*i + gjf) 
(hl + kl + h 2 k 2 + ^ 

) 


(26) 


(3) Twin Detection XRD for Cubic Crystalline Alloys on 
Trigonal/Hexagonal Substrate: 

Using the above equations (24), (25), and (26), the twin 
60 detection XRD method can be completed for all general cubic 
crystalline alloys on any trigonal or hexagonal substrate. 

As an example, calculated XRD peak positions for selected 
cubic crystalline alloy compositions on selected substrates 
65 are shown in the following tables. The actual compositions 
are not limited to these tables, but can be extended according 
to the equations from (3) to (20). 
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TABLE 1 
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TABLE 1 -continued 


Selected list of cubic ciystalline alloy compositions on various Selected list of cubic crystalline alloy compositions on various 

, . . , , .... . ^ ^ _ substrates in the lattice-match condition in Type A and Type B. 

substrates m the lattice-match condition in Type A and Type B. 1 1 


Epitaxial Layer 


Si Ge C Sn 




Alloy/Substrate- 

Lattice 





Type (A or B) 

C- of 

Substrate 

Lattice 


5.431 

5.658 

3.57 6.489 

pseudo 

Name 

C. 


Alloy Composition 

alloy 

Type A 

Type B 

Sapphire 

Substrate 







SiGe/ 

0.72 

0.28 


5.49 



Sapphire- A 
SiSn/ 

0.14 


0.86 

6.34 

Sapphire 

4.758 

Sapphire-B 

GeSn/ 


0.18 

0.82 

6.34 

5.49 

6.34 

Sapphire-B 

CSn/ 



0.05 0.95 

6.34 



Sapphire-B 

SiGeCSn/ 

Sapphire-B 

GaN 

Substrate 

0.03 

0.11 

0.01 0.85 

6.34 



SiC/GaN-A 

0.06 


0.94 

3.68 



GeC/GaN-A 


0.055 

0.945 

3.68 

GaN 

3.189 

SiC/GaN-B 

0.365 


0.635 

4.25 

3.68 

4.25 

GeC/GaN-B 


0.325 

0.675 

4.25 



SnC/GaN-B 



0.768 0.232 

4.25 



SiGeCSn/ 

0.11 

0.09 

0.7 0.1 

4.25 



GaN-B 







AIN/ 







Substrate 







SiC/AlN-A 

0.01 


0.99 

3.59 



SiC/AlN-B 

0.31 


0.69 

4.15 

AIN 

3.11 

GeC/AlN-B 


0.28 

0.72 

4.15 

3.59 

4.15 

CSn/AlN-B 



0.8 0.2 

4.15 




Epitaxial Layer 




Si 

Ge 

C 

Sn 







Alloy/Substrate' 


Lattice 






Type (A or B) 


C- of 

Substrate 

Lattice 

10 


5.431 

5.658 

3.57 

6.489 

pseudo 

Name 

C. 



Alloy Composition 

alloy 

Type A 

Type B 


SiGeC/ 

A1N-B 

0.2 

0.1 

0.7 


4.15 




SiGeCSn/ 

0.05 

0.05 

0.77 

0.13 

4.15 



15 

A1N-B 

InN 

Substrate 









SiC/InN-A 

0.28 


0.72 


4.09 




GeC/InN-A 


0.25 

0.75 


4.09 

InN 

3.54 

20 

CSn/InN-A 



0.822 

0.178 

4.09 

4.09 

4.72 

SiGeCSn/ 

InN-A 

0.05 

0.05 

0.79 

0.11 

4.09 





SiC/InN-B 

0.62 


0.38 


4.72 




GeC/InN-B 


0.55 

0.45 


4.72 




CSn/InN-B 



0.605 

0.395 

4.72 




SiGeCSn/ 

0.25 

0.05 

0.5 

0.2 

4.72 



25 

InN-B 

ZnO 

Substrate 









SiC/ZnO-A 

0.095 


0.905 


3.75 




GeC/ZnO-A 


0.084 

0.916 


3.75 

ZnO 

3.25 

30 

SiC/ZnO-B 

0.41 


0.59 


4.33 

3.75 

4.33 


CSn/ZnO-B 



0.74 

0.26 

4.33 




SiGeCSn/ 

ZnO-B 

Quartz 

Substrate 

0.195 

0.05 

0.655 

0.1 

4.33 



35 

Ge(Sn)/ 

Quartz-A 


0.99 


0.01 

5.67 




SiSn/ 

0.77 



0.23 

5.67 

Quartz 

4.91 


Quartz-A 

CSn/ 



0.28 

0.72 

5.67 

5.67 

6.55 

40 

Quartz-A 

SiGeCSn/ 

Quartz-A 

0.2 

0.2 

0.15 

0.45 

5.67 




TABLE 2 


Selected list of cubic crystalline alloy compositions on various substrates in the fractional 
lattice-match condition in 2 A of type A and V 2 of type B. 

Epitaxial Laver 



Si 

Ge 

C Sn 

Lattice 

Substrate 

Lattice 


Allov/Substrate- Typ e 

C. of 

Name 

C. 


5.431 

5.658 

3.57 6.489 

pseudo 

2 /3 Of 

V 2 of 


Alloy Composition 

alloy 

Type A 

Type B 

Sapphire Substrate 







SiC/Sapphire- 2 /3 of A 

0.05 


0.95 

3.66 

Sapphire 

4.758 

GeC/Sapphire- 2 /3 of A 
Lithium Niobate 
(LiNb03) Substrate 


0.045 

0.955 

3.66 

3.66 

3.17 

SiC/LiNb03- 2 / 3 of A 

0.21 


0.79 

3.96 

LiNb03 

5.148 

GeC/LiNb03- 2 / 3 of A 
Lithium Tantalum Nitride 
Oxide (LTNO) Substrate 


0.185 

0.815 

3.96 

3.96 

3.43 

SiC/LTNO - 2 / 3 of A 

0.555 


0.445 

4.60 

Lil6N801Ta2 

5.982 

SiC/LTNO - Vi of B 

0.225 


0.775 

3.99 

4.60 

3.99 
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TABLE 2 -continued 


Selected list of cubic crystalline alloy compositions on various substrates in the fractional 
lattice-match condition in 2 / 3 of type A and Vi of type B. 


Epitaxial Layer 





Si Ge C Sn 

Allov/Substrate-Type 

Lattice 
C. of 

Substrate 

Name 

Lattice 

C. 


5.431 5.658 3.57 6.489 

Alloy Composition 

pseudo 

alloy 

2 / 3 of 

Type A 

V 2 of 
Type B 

GeC/LTHO - Vi of B 
SiGeC/LTNO - Vi of B 
LiCuTa03 Substrate 
(Li = 0.667, Cu = 0.333) 

0.2 0.8 
0.1 0.11 0.79 

3.99 

3.99 



SiC/LiCu03 - 2 / 3 of A 
GeC/LiCu03 - 2 / 3 of A 

0.23 0.77 

0.205 0.795 

4.00 

4.00 

LiCuTa03 

4.00 

5.2004 

3.47 


The Bragg angles of X-ray diffraction were calculated and 
are listed in the following tables for select lattice matching 
conditions on various substrates. Table 3 and Table 5 were 
calculated with copper Cu K-cq X-ray source and Table 4 and 
Table 6 were calculated with X-ray from cobalt, with 
weighted average of Co K-a x and K-a 2 . 


TABLE 4 -continued 


Bragg angles (°) for selected XRD peaks of lattice matched alloy on 
various substrates (X-ray source: Co K-a [weighted average]). 

Co K-a (average) 

X-ray Wavelength: 1 .79026 Angstrom 
25 


TABLE 3 


1st order Bragg angle, Theta, for various peaks (in degrees) 


Bragg angles (°) for selected XRD peaks of lattice matched 
alloy on various substrates (X-ray source: Cu K-a^. 

Cu K-al 

X-ray Wavelength: 1.54056 Angstrom 

1st order Bragg angle. Theta, for various peaks (in degrees) 
L.C. of 



alloy 

{HI} 

{220} 

{113} 

{224} 

{004} 

Sapphire 

5.49 

14.065 

23.381 

27.732 

43.421 

34.140 

Type A 
Sapphire 

6.34 

12.148 

20.099 

23.763 

36.527 

29.077 

Type B 
GaN Type A 

3.68 

21.257 

36.301 

43.965 

N/A 

56.852 

GaN Type B 

4.25 

18.296 

30.839 

36.950 

62.611 

46.466 

AIN Type A 

3.59 

21.816 

37.364 

45.367 

N/A 

59.121 

AIN Type A 

4.15 

18.753 

31.667 

37.995 

65.408 

47.940 

InN Type A 

4.09 

19.038 

32.187 

38.655 

67.315 

48.880 

InN Type B 

4.72 

16.419 

27.489 

32.769 

53.081 

40.752 

ZnO Type A 

3.75 

20.841 

35.520 

42.942 

N/A 

55.248 

ZnO Type B 

4.33 

17.946 

30.209 

36.157 

60.633 

45.363 

Quartz 

5.67 

13.609 

22.597 

26.780 

41.723 

32.916 


Type A 


TABLE 4 


Bragg angles (°) for selected XRD peaks of lattice matched alloy on 
various substrates (X-ray source: Co K-a [weighted average]). 

Co K-a (average) 

X-ray Wavelength: 1.79026 Angstrom 

1st order Bragg angle. Theta, for various peaks (in degrees) 
L C al- 



loy 

{in} 

{220} 

{113} 

{224} 

{004} 

Sapphire 

5.49 

16.404 

27.462 

32.736 

53.012 

40.707 

Type A 
Sapphire 

6.34 

14.155 

23.537 

27.922 

43.763 

34.385 

Type B 
GaN Type A 

3.68 

24.917 

43.471 

53.779 

N/A 

76.647 

GaN Type B 

4.25 

21.395 

36.564 

44.310 

N/A 

57.402 

AIN Type A 

3.59 

25.586 

44.849 

55.788 

N/A 

85.835 

AIN Type A 

4.15 

21.937 

37.595 

45.674 

N/A 

59.630 

InN Type A 

4.09 

22.276 

38.245 

46.541 

N/A 

61.096 


30 


35 


40 


45 


50 


55 


60 


65 


L C_al- 

loy {111} {220} {113} {224} {004} 


InN Type B 

4.72 

19.176 

ZnO Type A 

3.75 

24.421 

ZnO Type B 

4.33 

20.981 

Quartz 

5.67 

15.869 

Type A 




32.439 

38.975 

68.291 

49.340 

42.466 

52.343 

N/A 

72.709 

35.783 

43.286 

N/A 

55.783 

26.521 

31.574 

50.661 

39.160 


TABLE 5 


Bragg angles (°) for selected XRD peaks of alloys 
of fractional lattice matching condition on various 
substrates (X-ray source: Cu K-a t ). 

1st order Bragg angle. Theta, for various peaks (in degree) 


L.C._al- 

loy {111} 


Sapphire - 

2 / 3 of 

3.66 

21.378 

Type A 
Sapphire - 
Vi of 

3.17 

24.890 

Type B 
LiNb03 - 

3.96 

19.689 

2 / 3 of 



Type A 
LiNb03 - 

3.43 

22.890 

V 2 of 
Type B 
LTNO - 

2 / 3 of 

4.6 

16.860 

Type A 
LTNO - 

3.99 

19.535 

V 2 of 
Type B 
LiCuTa03 - 

2 / 3 of 

4.00 

19.484 

Type A 
LiCuTa03 - 
V 2 of 
Type B 

3.47 

22.612 


{220} 

{113} 

{224} 

{004} 

36.532 

44.268 

N/A 

57.335 

43.415 

53.698 

N/A 

76.400 

33.379 

40.176 

N/A 

51.083 

39.434 

48.144 

N/A 

63.934 

28.270 

33.737 

N/A 

42.052 

33.095 

39.813 

71.043 

50.553 

33.002 

39.694 

70.630 

50.379 

38.893 

47.412 

N/A 

62.614 
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TABLE 6 


Bragg angles (°) for selected XRD peaks of alloys of 
fractional lattice matching condition on various substrates 
(X-ray source: Co K-a [weighted average!) 

1st order Bragg angle. Theta, for various peaks (in degree) 

L.C._al- 



loy 

{ 111 } 

{220} 

{113} 

{224} 

{004} 

Sapphire - 

2 /3 Of 

3.66 

25.063 

43.769 

54.208 

N/A 

78.038 

Type A 
Sapphire - 
Viof 

3.17 

29.281 

53.004 

69.477 

N/A 

N/A 

Type B 
LiNb03 - 

2 /3 Of 

3.96 

23.049 

39.743 

48.564 

N/A 

64.712 

Type A 
LiNb03 - 
Viof 

3.43 

26.873 

47.573 

59.945 

N/A 

N/A 

Type B 
LTNO - 

4.6 

19.697 

33.394 

40.195 

72.423 

51.112 

2/3 Of 

Type A 
LTNO - 
Viof 

3.99 

22.866 

39.386 

48.079 

N/A 

63.815 

Type B 
LiCuTa03 - 

2/3 of 

4 

22.805 

39.268 

47.919 

N/A 

63.525 

Type A 
LiCuTa03 - 

3.47 

26.539 

46.855 

58.822 

N/A 

N/A 


Vi of 

Type B 


The inter-planar angle between {111} plane and a few 
other planes are listed in the following Table 7. 


TABLE 7 


Inter-planar angle between f 1 1 1| plane and various other planes. 


Inter-Planar 

Planes 

Angle (°) 

{111}/ {224} 

19.471 

{111 }/{ 113} 

29.496 

{111}/ {220} 

35.264 

{111}/ {004} 

54.736 


These calculated numbers can be applied to the technology 
of copending application Ser. No. 12/254,150, entitled 
“METHOD OF GENERATING X-RAY DIFFRACTION 
DATA FOR INTEGRAL DETECTION OF TWIN 
DEFECTS IN SUPER-HETERO-EPITAXIAL MATERI- 
ALS,” and Ser. No. 12/288,380, entitled “X-RAY DIFFRAC- 
TION WAFER MAPPING METHOD FOR RHOMBOHE- 
DRAL SUPER-HETERO-EPITAXY, ” which are both 
incorporated herein by reference in their entirety. 

Therefore, the twin detection XRD can be expanded to all 
group IV alloys, as well as group III-V and group II- VI 
crystalline materials in zinc -blende structure on any trigonal 
or hexagonal substrate, with the Bragg angles that are calcu- 
lated in Tables 2 and 3 according to equation (19) and the 
inter-planar angles that are calculated in the Table 4 according 
to the equation (20). With appropriate [h,j,k] values, the same 
techniques could also be applied to other BCC and FCC 
materials. 

The equations and tabulated numbers above can be used to 
apply the methods of accompanying application Ser. No. 
1 2/254,1 50, entitled “METHOD OF GENERATING X-RAY 
DIFFRACTION DATA FOR INTEGRAL DETECTION OF 
TWIN DEFECTS IN SUPER-HETERO-EPITAXIAL 
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MATERIALS,” and Ser. No. 12/288,380, entitled “X-RAY 
DIFFRACTION WAFER MAPPING METHOD FOR 
RHOMBOHEDRAL SUPER-HETERO-EPITAXY” which 
are both incorporated herein by reference in their entirety, for 
5 XRD twin detection to all group IV semiconductor alloys 
made from any of Si, Ge, C, and Sn elements, and group III-V 
and group II -VI crystalline materials in zinc -blende structure, 
on any trigonal/hexagonal substrate. 

io ADVANTAGES OF THE TECHNOLOGY 

The rhombohedral-trigonal (hexagonal) lattice matching 
technology described herein enables the development of new 
cubic crystalline semiconductor alloys substantially relieved 
15 from misfit dislocation defects on various substrates. One 
example is shown in FIG. 10 for today’s fastest transistor, 
Hetero Bipolar Transistor (HET), which is made with SiGe 
alloy with small Ge content. The electron mobility of the thin 
Base region of HBT is much faster than that of conventional 
20 Si transistors because the conduction band of the gradient 
SiGe alloy layer in the Base region has a built-in electric field 
due to the band-gap energy difference between the Si layer 
and the SiGe layer, which is shown as “Slope” in FIG. 10. In 
spite of these advantages, it has been extremely difficult to 
25 make a defect- free SiGe alloy with high Ge content because 
the lattice constant of Ge is larger than that of Si, as shown at 
the top of FIG. 10. This difference resulted in a critical thick- 
ness, which is a theoretical limit that epitaxial layer can be 
grown without dislocation defects. 

30 The present invention makes it possible to incorporate high 

content of Ge and Tin (Sn) atoms into Si free from the prior 
critical thickness limit because it is lattice-matched, in order 
to increase the built-in electric field for even faster SiGeCSn 
HBT devices and more efficient SiGeCSn TE devices. Note 
35 the greater slope in the graph for SiGeCSn on the lower 
right-hand side of FIG. 10. Because faster mobility of elec- 
trons and holes are very important to increase ZT factor of TE 
materials, this technology will be of great use in developing 
faster semiconductors. 

40 Thus, it is apparent that the methods and materials 
described in the present disclosure satisfy the objects of the 
invention set forth above. 

The present invention is further discussed in Y. Park, G. C. 
King, S. H. Choi, Rhombohedral epitaxy of cubic SiGe on 
45 trigonal c -plane sapphire, Journal of Crystal Growth 310 
(2008) 2724-2731, herein incorporated by reference in its 
entirety. 

Although the present invention has been described in 
detail, it should be understood that various changes, substi- 
50 tutions, and alterations may be readily ascertainable by those 
skilled in the art and may be made herein without departing 
from the spirit and scope of the present invention as defined 
by the following claims. 

55 What is claimed as new and desired to be secured by 
Letters Patent of the United States is: 

1. A method of forming a hetero-epitaxial crystal structure 
of a group IV, group III-V or group II-VI alloy on the basal 
(0001) plane of a trigonal crystal or hexagonal crystal sub- 
60 strate so as to minimize crystalline dislocation defects 
between the hetero-epitaxial crystal structure of the group IV, 
group III-V or group II-VI alloy and the trigonal crystal or 
hexagonal crystal substrate, said method comprising: 

(a) if said alloy is a type in which the trigonal arrangement 
65 of atomic sites of the {111} plane of said alloy crystal 
directly meets the hexagonal or trigonal coincidence 
lattice of atomic sites in the basal plane of said substrate 
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without any rotation, calculating a target lattice constant 
(L target) as 2A/3*basal plane lattice constant of the sub- 
strate (a m/Ae J; 

(b) if said alloy is of a type in which the trigonal arrange- 
ment of atomics sites of the {111} plane of said alloy 5 
crystal meets the hexagonal or trigonal coincidence lat- 
tice of atomic sites in the basal plane of said substrate 
with a 30° rotation, calculating E target as 4/3*basal plane 
lattice constant of the substrate (a tri/hex )', 

(c) if said alloy is of a type in which the trigonal arrange- 10 
ment of atomic sites of the {111} plane of said alloy 
crystal directly meets the inner vertices of the hexagonal 

or trigonal coincidence lattice of atomic sites in the basal 
plane of said substrate without any rotation, calculating 
a target lattice constant (E target ) as 4/ (3 *V3)*basal plane 15 
lattice constant of the substrate (a tri/hex )', 

(d) if said alloy to be of a type in which the trigonal 
arrangement of atomic sites of the {111} plane of said 
alloy crystal meets the inner vertices of the hexagonal or 
trigonal coincidence lattice of atomic sites in the basal 20 
plane of said substrate with a 30° rotation, calculating 

L target as %*basal plane lattice constant of the substrate 

fotri/hex) "> 

(e) determining the respective proportions of the group IV, 

group III-V or group II-VI elements in said alloy by 25 
approximating E target as a linear function of the indi- 
vidual lattice constants for the constituent group IV, 
group III-V or group II-VI elements weighted in accor- 
dance with said proportions and a Bowing Parameter, 
and solving for said proportions; and 30 

(f) growing the group IV, group III-V or group II-VI alloy 
on the basal (0001) plane of the trigonal crystal or hex- 
agonal crystal substrate in the determined respective 
proportions to thereby form the hetero -epitaxial crystal 
structure on the trigonal crystal or hexagonal crystal 
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substrate so as to minimize crystalline dislocation 
defects between the hetero-epitaxial crystal structure of 
the group IV, group III-V or group II-VI alloy and the 
trigonal crystal or hexagonal crystal substrate. 

2. The method of claim 1, wherein the Bowing Parameter is 
treated as zero, thereby reducing step c to the step of solving 
a system of linear equations. 

3. The method of claim 2, wherein said alloy is a binary 
alloy in which the two group IV, group III-V or group II-VI 
components have relative proportions of x and 1 -x and lattice 
constant L A and L B , and, after E target has been determined in 
accordance with steps (a) or (b) in claim 1, x is determined by 
solving the equation x*L^+(l-x)*L^=L to ^ er 

4. The method of claim 2, wherein said alloy is a ternary 
alloy in which the three group IV, group III-V or group II-VI 
components have relative proportions of x, y and 1-x-y, 
lattice constants L a ,L 3 and L c , and band-gap energies E A , E B 
and E c , and, after E target has been determined in accordance 
with steps (a) or (b) in claim 1, x and y are determined by 
solving the equations x* +y * Lg+ ( 1 -x-y) * E c =L target and 
x*E A +y*E B +(l-x-y)*E c =E target and the solution is chosen 
for which x and y are each between 0 and 1 . 

5. The method of claim 2, wherein said alloy is a quaternary 
alloy in which the four group IV, Group III-V or group II-VI 
components have relative proportions of x, y, z and 1 -x-y-z, 
lattice constants L A , L B , L c and L D , band-gap energies E A , E B , 
E c and E D , and linear thermal expansion coefficients cx 4 , a B , 
a c , and a D , and, after E target has been determined in accor- 
dance with steps (a) or (b) in claim 1, x, y and z are determined 
by solving the equations x*L^+y*L 5 +z*L c +(l-x-y-z) 
*Lr.=L MJSep x*E A +y*E B +z*E c +(l-x-y-z)*E n =E target , and 
x*a A +y*a B + z *a.c + (^- x -y- z )* a n =0i target> an d the solution 
is chosen for which x, y and z are each between 0 and 1 . 



